Bose-Einstein condensation of polaritons in graphene in a high magnetic field 
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The Bose-Einstein condensation (BEC) of magnetoexcitonic polaritons in a graphene layer em- 
bedded in a optical microcavity in a high magnetic field B is predicted. The essential property of this 
system (in contrast, e.g., to a quantum well embedded in a cavity) is stronger influence of magnetic 
field and weaker influence of disorder. A two-dimensional (2D) magnetoexcitonic polaritons gas is 
considered in a planar harmonic electric field potential applied to excitons or a parabolic shape of 
the optical cavity causing the trapping of microcavity photons. It is shown that the effective polari- 
ton mass M e g increases with magnetic field as B 1 ^ 2 . The BEC critical temperature Ti°' decreases 
as B 1,/4 and increases with the spring constant of the parabolic trap. The Rabi splitting related 
to the creation of a magnetoexciton in a high magnetic field in graphene is obtained. 

PACS numbers: 71.36.+C, 03.75.Hh, 73.20.Mf, 73.21.Fg 

In the past decade, Bose coherent effects of 2D excitonic polaritons in a quantum well embedded in a semiconductor 
microcavity have been the subject of theoretical and experimental studies To obtain polaritons, two mirrors placed 
opposite each other form a microcavity, and quantum wells are embedded within the cavity at the antinodes of the 
(_J . confined optical mode. The resonant exciton-photon interaction results in the Rabi splitting of the excitation spectrum. 
Two polariton branches appear in the spectrum due to the resonant exciton-photon coupling. The lower polariton 
(LP) branch of the spectrum has a minimum at zero momentum. Recently, the polaritons in a harmonic potential 
trap have been studied experimentally in a GaAs/AlAs quantum well embedded in a GaAs/AlGaAs microcavity [H- 
In this trap, the exciton energy is shifted using a stress-induced band-gap. In this system, evidence for the BEC of 
polaritons in a quantum well has been observed Q ■ The theory of the BEC and superfluidity of exciton polaritons in 
I a quantum well (QW) without magnetic field in a parabolic trap has been developed in Ref. [3| . 

A novel type of 2D electron system was experimentally obtained in graphene, which is a 2D honeycomb lattice of 
the carbon atoms that form the basic planar structure in graphite 0, Q . Graphene is a gapless semiconductor with 
massless electrons and holes which have been described as Dirac-fermions. The energy spectrum and the wavefunctions 
t^J- ' of electron-hole pairs in a magnetic field, in graphene have been calculated in Ref. [1, 0). In high magnetic field 
£C) , electron-hole pairs in graphene form two-dimensional magnetoexcitons. The BEC and superfluidity of spatially- 
indirect magnetoexcitons with spatially separated electrons and holes have been studied in graphene bilayer Q and 
graphene superlattice @. The electron-hole pair condensation in the graphene-based bilayers have been studied 
■ in [II, E2, E2, EH . The similar effects in the system of spatially- indirect excitons (or electron-hole pairs) in coupled 
quantum wells (CQWs), with and without a magnetic field were studied in Refs. [3, El El, El 03 EH- The 
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ON , experimental and theoretical interest to study these systems is particularly due to the possibility of the BEC and 
superfluidity of indirect excitons, which can manifest in the CQW as persistent electrical currents in each well and 
also through coherent optical properties and Josephson phenomena 0, El, EH E3, El, EH • The great experimental 
success was achieved now in this field [H, [2l|, [U HI]- Besides, the essential experimental progress was achieved 
in experimental studies of exciton polaritons in the system of a QW embedded in optical microcavity [2J, |25|, |26| . 
However, while the exciton polaritons have been studied in a QW, the formation of the polaritons in graphene in high 
magnetic field have not yet been considered. Moreover, the polaritons formed as superposition of magnetoexcitons 
and cavity photons in magnetic field have not yet been studied. It is interesting to study a 2D system of polaritons 
in graphene embedded in a microcavity from the point of view of the existence of the BEC within it. The BEC of 
polaritons in high magnetic field in graphene embedded in a microcavity is perspective, since the random field in 
graphene is weaker than in a QW, particularly, because in a QW the random field is generated due to the fluctuations 
of the width of a QW. Let us mention that if the interaction of bosons with the random field is stronger, the BEC 
critical temperature is lower [27] . Moreover, the BEC of polaritons in graphene embedded in a microcavity can exist 
at much lower magnetic field than in a QW, because the distance between electron Landau levels in graphene is much 
higher than in a QW at the same magnetic field, and, therefore, the lower magnetic field is required in graphene than 
in a QW to neglect the electron transitions between the Landau levels. The purpose of this Letter is to point out the 
existence of the BEC of the magnetoexciton polaritons in a graphene layer embedded in a semiconductor microcavity 
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in a strong magnetic field and to discuss the conditions of its realization. 

When an undoped electron system in graphene in a magnetic field without an external electric field is in the ground 
state, half of the zeroth Landau level is filled with electrons, all Landau levels above the zeroth one are empty, and all 
levels below the zeroth one are filled with electrons. We suggest using the gate voltage to control the chemical potential 
in graphene by two ways: to shift it above the zeroth level so that it is between the zeroth and first Landau levels (the 
first case) or to shift the chemical potential below the zeroth level so that it is between the first negative and zeroth 
Landau levels (the second case). In both cases, all Landau levels below the chemical potential are completely filled 
and all Landau levels above the chemical potential are completely empty. Based on the selection rules for optical 
transitions between the Landau levels in single-layer graphene [28|], in the first case, there are allowed transitions 
between the zeroth and the first Landau levels, while in the second case there are allowed transitions between the 
first negative and zeroth Landau levels. Correspondingly, we consider magnetoexcitons formed in graphene by the 
electron on the first Landau level and the hole on the zeroth Landau level (the first case) or the electron on the zeroth 
Landau level and the hole on the Landau level —1 (the second case). Note that by appropriate gate potential we can 
also use any other neighboring Landau levels n and n + 1. 

For the relatively high dielectric constant of the microcavity, e ^ e 2 /(hvF) ~ 2 (vf = V3at/(2H) is the Fermi 
velocity of electrons in graphene, where a = 2.566 A is a lattice constant and t ~ 2.71 eV is the overlap integral 
between the nearest carbon atoms (29j ) the magnetoexciton energy in graphene can be calculated by applying 
perturbation theory with respect to the strength of the Coulomb electron-hole attraction analogously as it was done 
in [3p| for 2D quantum wells in a high magnetic field with non-zero electron m e and hole ro/j masses. This approach 
allows us to obtain the spectrum of an isolated magnetoexciton with the electron on the Landau level 1 and the hole 
on the Landau level in a single graphene layer, and it will be exactly the same as for the magnetoexciton with 
the electron on the Landau level and the hole on the Landau level —1. The characteristic Coulomb electron- hole 
attraction for the single graphene layer is e 2 /(ere), where e is the dielectric constant of the environment around 
graphene, tb = -J ch/(eB) denotes the magnetic length of the magnetoexciton in the magnetic field B, and c is 
the speed of light. The energy difference between the first and zeroth Landau levels in graphene is Iivf/tb- For 
graphene, the perturbative approach with respect to the strength of the Coulomb electron-hole attraction is valid 
when e 2 /(ere) <C YivfItb [30|- This condition can be fulfilled at all magnetic fields B if the dielectric constant of the 
surrounding media satisfies the condition e 2 /(ehvp) <C 1. Therefore, we claim that the energy difference between the 
first and zeroth Landau levels is always greater than the characteristic Coulomb attraction between the electron and 
the hole in the single graphene layer at any B if e ^> e 2 /(Hvf) ~ 2. Thus, applying perturbation theory with respect 
to weak Coulomb electron-hole attraction in graphene embedded in the GaAs microcavity (e = 12.9) is more accurate 
than for graphene embedded in the SiC>2 microcavity (e = 4.5). However, the magnetoexcitons in graphene exist in 
high magnetic field. Therefore, we restrict ourselves by consideration of high magnetic fields. 

Polaritons are linear superpositions of excitons and photons. In high magnetic fields, when magnetoexcitons may 
exist, the polaritons become linear superpositions of magnetoexcitons and photons. Let us define the superpositions 
of magnetoexcitons and photons as magnetopolaritons. It is obvious that magnetopolaritons in graphene are two- 
dimensional, since graphene is a two-dimensional structure. The Hamiltonian of magnetopolaritons in the strong 
magnetic field is given by 

P-tot — Hrnex P-ph ~t~ Prnex—ph j (1) 

where H mex is a magnetoexcitonic Hamiltonian, H p h is a photonic Hamiltonian, and H exc ^ p h is the Hamiltonian of 
magnetoexciton-photon interaction. 

Let us analyze each term of the Hamiltonian for magnetopolaritons JT]). The effective Hamiltonian and the energy 
dispersion for magnetoexcitons in graphene layers in a high magnetic field derived in Ref. Q are given by 

(P)bpb P , e mex (P) = E band - £^> + e (P) , (2) 

p 

where bp and bp are magnetoexcitonic creation and annihilation operators obeying the Bose commutation relations. 
In Eq. @, Eb anc i — e[°q — \f2hvF Itb is the band gap energy, which is the difference between the Landau levels 1 and 

in graphene. is the binding energy of a 2D magnetoexciton with the electron in the Landau level 1 and the hole 
on the Landau level in a single graphene layer, and £o{P) = P 2 /(2itlb), where raj is the effective magnetic mass 
of a 2D magnetoexciton with the electron on the Landau level 1 and hole on the Landau level in a single graphene 
layer. The binding energy £^ and effective magnetic mass tub of a magnetoexciton in graphene obtained using the 
first order perturbation respect to the electron-hole Coulomb attraction similarly to the case of a single quantum well 
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[301 ] are given by 



Ab ) * e 2 2 7 / 2 eft 2 



E B =\ o ' mB = —i=-2 • ( 3 ) 



It can be shown that the interaction between two direct 2D magnetoexcitons in graphene with the electron on the 
Landau level 1 and the hole on the Landau level can be neglected in a strong magnetic held, in analogy to what is 
described in Ref . [30l ] for 2D magnetoexcitons in a quantum well. The dipole moment of each exciton in a magnetic field 
is d 12 = ep = t\ [B x Pi. 2] /B [30], where Pi and P 2 are the magnetic momenta of each exciton and Pi, P2 <C 
The magnetoexcitons are located at a distance i? 3> tb from each other. The corresponding contribution to the energy 
of their dipole-dipole interaction is - £^ (r B /R) 3 PiP 2 r 2 B /e ~ (r B /Rf P 1 P 2 /(eM ) < e 2 r|/(ei? 3 ). Inputting the 
radius of the magnetoexciton in graphene ro,i ~ Tb 0, we obtain that the van der Waals attraction of the exciton 
at zero momenta is proportional to ~ (r 0j i/i?) 6 ~ (r^/i?) 6 . Therefore, in the limit of a strong magnetic field for 
a dilute system tb <C R, both the dipole-dipole interaction and the van der Waals attraction vanish, and the 2D 
magnetoexcitons in graphene form an ideal Bose gas analogously to the 2D magnetoexcitons in a quantum well given 
in Ref. [30(. Thus, the Hamiltonian |T]) does not include the term corresponding to the interaction between two direct 
magnetoexcitons in a single graphene layer. So in high magnetic field there is the BEC of the ideal magnetoexcitonic 
gas in graphene. 

The Hamiltonian and the energy spectrum of non-interacting photons in a semiconductor microcavity are given by 



H ph = J2 s ph{P)a±>a P , e ph {P) = (c/n^P 2 + h^ 2 L c 2 , (4) 
p 

where a P and ap are photonic creation and annihilation Bose operators. In Eq. Lc is the length of the cavity, 
n = y / £c is the effective refractive index and ec is the dielectric constant of the cavity. We assume that the length of 
the microcavity has the following form: 

/ „n hire , , 

Lc{B) = ~Tf — m ' (5) 

corresponding to the resonance of the photonic and magnetoexcitonic branches at P = 0, i.e. e mex (0) — £ p /j(0). 
As it follows from the energy spectra in and (j4]), and Eqs. ((3]) and {S]), the length of the microcavity, corre- 
sponding to a magnetoexciton-photon resonance, decreases with the increment of the mag netic field as B^ 1 / 2 . The 
resonance between magnetoexcitons and cavity photonic modes can be achieved either by controlling the spectrum of 
magnetoexcitons e e x{P) by changing magnetic field B or by choosing the appropriate length of the microcavity Lc- 
The Hamiltonian of the magnetoexciton-photon coupling has the form (see Refs. [32T l33l l34j|): 

Hmex-ph = HCIr apbp + h.c. , (6) 



bp = Xpp-p — Cpu-p, ap = Cpp-p + Xpup , (7) 

where the magnetoexciton-photon coupling energy represented by the Rabi constant Mlp is obtained below, pp and 
■up are lower and upper magnetopolariton Bose operators, respectively, j Xp and Cp are coefficients of the unitary 
Bogoliubov transformation [32|, [34|, and the energy spectra of the lower/upper magnetopolaritons are 

fm Sph(P) + £mex(P) 

£lp/up{P) = 7, 



T ^(£ P h(P) - e mex {P)f + 4|M7 fl | 2 . (8) 

Eq. ([5]) implies a splitting of 2ft£lp between the upper and lower states of polaritons at P = 0, which is known as the 
Rabi splitting. 

Substituting into Hamiltonians (O, (UJ) and ^ and diagonalizing the Hamiltonian |T]), finally we obtain 

H tot = £ Lp{P)PpPP + £up(P)upu P , (9) 
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where pp, pp, up, up are the Bose creation and annihilation operators for the lower and upper magnetopolaritons, 
respectively. The Hamiltonian Eq. ^ describes magnetopolaritons in a single graphene layer in a high magnetic field. 
Our particular interest is the lower energy magnetopolaritons which produce the BEC. The lower palaritons have the 
lowest energy within a single graphene layer. 

Similarly to the Bose atoms in a trap in the case of a slowly varying external potential (35j . we can make the 
quasiclassical approximation, assuming that the effective magnetoexciton mass does not depend on a characteristic 
size I of the trap and it is a constant within the trap. This quasiclassical approximation is valid if P 3> h/l. The 
harmonic trap is formed by the two-dimensional planar potential in the plane of graphene. The potential trap can be 
produced in two different ways. In case 1, the potential trap can be produced by applying an external inhomogencous 
electric field. The spatial dependence of the external field potential V(r) is caused by shifting of magnetoexciton 
energy by applying an external inhomogeneous electric field. The photonic states in the cavity are assumed to be 
unaffected by this electric field. In this case the band energy Et an d is replaced by Eb an d(r) — -Eband(O) + V(r). Near 
the minimum of the magnetoexciton energy, V(r) can be approximated by the planar harmonic potential jr 2 /2, 
where 7 is the spring constant. Note that a high magnetic field does not change the trapping potential in the effective 
Hamiltonian [361 ] . In case 2, the trapping of magnetopolaritons is caused by the inhomogeneous shape of the cavity 
when the length of the cavity is determined by Eq. (0 with the term jr 2 /2 added to Eb an d — £g\ where r is the 
distance between the photon and the center of the trap. In case 2, the 7 is the curvature characterizing the shape 
of the cavity. In case 1, for the slowly changing confining potential V(r) = 7r 2 /2, the magnetoexciton spectrum is 

given in the effective mass approximation by Eq. @ provided we added the term r yr 2 /2 added to Eb an d — £^ in 
the r.h.s., where r is now the distance between the center of mass of the magnetoexciton and the center of the trap. 
The Hamiltonian and the energy spectrum of the photons in this case are shown by Eq. ^ , and the length of the 
microcavity is given by Eq. ([5]). 

In case 2, for the photonic spectrum in the effective mass approximation is given by substituting the slowly changing 
shape of the length of cavity depending on the term 7r 2 /2 into Eq. ^ representing the spectrum of the cavity photons. 

1 /2 

This quasiclassical approximation is valid if P 3> H/l, where / = (H/ (rriBWo)) is the size of the magnetoexciton 
cloud in an ideal magnetoexciton gas and ujq = \J j/ms- The Hamiltonian and energy spectrum of magnetoexcitons 
in this case are given by 



At small momenta a = \/2{m B 1 + [c/n)Lc /hm)P 2 /\MIr\ <C 1 ( Lq = fmc/n [Eband — Eg] ) and weak confine- 



i -1 

l B T^/u^c/raijj / I'M'iil ^ J- V — ii'i'L./ ii, 1 ^band — o fl 

ment (3 = jr 2 /\hfln\ <C 1, the single-particle lower magnetopolariton spectrum obtained through the substitution of 
Eq. ([2]) into Eq. ([8]), in linear order with respect to the small parameters a and f3, is 

e (P) « ^HttLq 1 - \m R \ + lr 2 + \ (rn? + ^) P 2 . (10) 

Let us emphasize that the spectrum of non- interacting magnetopolaritons £q(P) at small momenta and weak confine- 
ment is given by Eq. (|10[) for both physical realizations of confinement: case 1 and case 2. 

If we measure the energy relative to the P = lower magnetopolariton energy (c/rijhTrL^, 1 — |M1r|, we obtain the 
resulting effective Hamiltonian for trapped magnetopolaritons in graphene in a magnetic field. At small momenta 

a <C 1 {Lc = hnc/n \ Eband ~ ^bj ) ano - wea k confinement /3< 1, this effective Hamiltonian is 

where the sum over P is carried out only over P ^> H/l (only in this case the quasiclassical approach used in Eq. (|2|) 
is valid), and the effective magnetic mass of a magnetopolariton is given by 

M .„ (B) = 2 ( mil + £M£>)-\ (12) 

According to Eqs. (fT2)) and ([3]), the effective magnetopolariton mass M e g increases with the increment of the magnetic 
field as B 1 ! 2 . Let us emphasize that the resulting effective Hamiltonian for magnetopolaritons in graphene in a 
magnetic field for the parabolic trap is given by Eq. (jlip for both physical realizations of confinement represented by 
case 1 and case 2. 

Neglecting anharmonic terms for the magnetoexciton-photon coupling, the Rabi splitting constant can be 
estimated quasiclassically as 

\Ml R \ 2 = ( f H int i) , Epho = I „ r I , H int = o ■ A= - — <7 • E ph0 , (13) 

/ \ eW I c icu 
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where a = (a x ,a y ), a x and a y are Pauli matrices, Hi nt is the Hamiltonian of the electron-photon interaction cor- 
responding to the electron in graphene described by Dirac dispersion, is the electric field corresponding to a 
single cavity photon, W is the volume of microcavity, u> is the photon frequency. The initial \i) electron state cor- 
responds to the completely filled Landau level and completely empty Landau level 1. The final |/) electron state 
corresponds to creation of one magnetoexciton with the electron on the Landau level 1 and the hole on the Landau 
level 0. The transition dipole moment corresponding to the process of creation of this magnetoexciton is given by 
di2 = ere/4. Let us note that in Eq. (fT3")) the energy of photon absorbed at the creation of the magnetoexciton is given 

by hw = E\ — £q = V^hvp/rs (we assume that £g <C E\ — Eq). Substituting the photon energy and the transition 
dipole moment from into Eq. (| 13[) , we obtain the Rabi splitting corresponding to the creation of a magnetoexciton with 

1 /2 

the electron on the Landau level 1 and the hole on the Landau level in graphene: KIr = e (irhvprB /{\/2eW)) 

Thus, the Rabi splitting in graphene is related to the creation of the magnetoexciton, which decreases when the 
magnetic field increases and is proportional to i? -1 / 4 . Therefore, the Rabi splitting in graphene can be controlled by 
the external magnetic field. It is easy to show that the Rabi splitting related to the creation of the magnetoexciton, the 
electron on the Landau level and the hole on the Landau level — 1 will be exactly the same as for the magnetoexciton 
with the electron on the Landau level 1 and the hole on the Landau level 0. Let us mention that dipole optical 
transitions from the Landau level —1 to the Landau level 0, as well as from the Landau level to the Landau level 1, 
are allowed by the selection rules for optical transitions in single-layer graphene [281 ]. 

As known, the exact solution for the ground state of 2D magnetoexcitons in high magnetic fields is the Bose 
condensate of non-interacting magnetoexcitons [30j |. This is valid for both semiconductor quantum well and graphene. 
Although Bose-Einstein condensation cannot take place in a 2D homogeneous ideal gas at non-zero temperature, as 
discussed in Ref. (37j . in a harmonic trap the BEC can occur in two dimensions below a critical temperature T®. In a 
harmonic trap at a temperature T below a critical temperature Xj? (T < T°), the number Nq(T, B) of non- interacting 
magnetopolaritons in the condensate is given in Ref. [37]. Applying the condition Nq = 0, and assuming that the 
magnetopolariton effective mass is given by Eq. (|12p. we obtain the BEC critical temperature for the ideal gas 
of magnetopolaritons in a single graphene layer in a magnetic field: 




T?\B) = JL I , - - 3 *7* , I , (14) 



1/2 

{gi e) gi e) +gi h) gi h) )M eB (B). 

where N is the total number of magnetopolaritons, gi and gi are the spin and graphene valley degeneracies 
for an electron and a hole, respectively, fee is the Boltzmann constant. At temperatures above Tc°\ the BEC of 
magnetopolaritons in a single graphene layer does not exist. Tc^/y/N as a function of magnetic field B and spring 
constant 7 is presented in Fig. [T] In our calculations, we used gs = gi^ = gs = g„ = 2. According to Eq. (|14|). 

the BEC critical temperature decreases with the magnetic field as _B~ 1//4 and increases with the spring constant 
as 7 1 / 2 . Note that we assume that the quality of the cavity is sufficiently high, so that the time of the relaxation to 
the Bose condensate quasiequilibrium state is smaller than the life time of the photons in the cavity. 

To conclude, we have studied the ideal gas of trapped cavity magnetopolaritons in a single graphene layer in a 
high magnetic field. The resonance between magnetoexcitons and cavity photonic modes can be achieved either by 
controlling the spectrum of magnetoexcitons e ex (P), by changing magnetic field B or by controlling the length of 
the microcavity Lq- We analyzed two possible physical realizations of the trapping potential: a harmonic electric 
field potential coupled to magnetoexcitons and a parabolic shape of the semiconductor cavity causing the trapping 
of microcavity photons. We conclude that both realizations of confinement result in the same effective Hamiltonian. 
It is shown that the effective magnetopolariton mass M e g increases with the magnetic field as B 1 / 2 . Meanwhile, 
the BEC critical temperature decreases as _B -1 / 4 and increases with the spring constant as 7 1 / 2 . The gas of 
magnetopolaritons in graphene in a high magnetic field can be treated as an ideal Bose gas since magnetoexciton- 
magnetoexciton interaction vanishes in the limit of a high magnetic field and a relatively high dielectric constant of 
the cavity e ^> 2. Let us mention that this condition for the high dielectric constant of the microcavity is valid only 
for graphene, and it is not valid for the quantum well. Besides, we have obtained the Rabi splitting related to the 
creation of a magnetoexciton in a high magnetic field in graphene which can be controlled by the external magnetic 
field B. 
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FIG. 1: The ratio of the BEC critical temperature to the square root of the total number of magnetopolaritons Tc^/y/N as 
a function of magnetic field B at different spring constants 7. We assume the environment around graphene is GaAs with 
e = 12.9. 
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